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LUNAR GRAVITATIONAL FIELD; A THIRTEENTH DEGREE 


AND ORDER SPHERICAL HARMONIC ESTIMATE 

By W. Thomas Blackshear, Edward F. Daniels, 
and S. G. Anderson 
Langley Research Center 

SUMMARY 

Lunar Orbiter tracking data have been analyzed to yield a thirteenth degree and 
order spherical harmonic approximation to the lunar gravitational potential function. 

The results of the analysis indicate an essentially homogeneous moon and demonstrate a 
correlation between lunar surface features and gravitational anomalies. In addition, a 
correlation exists between the type of data coverage and the general character of accel- 
eration anomalies calculated from the potential approximation. 

INTRODUCTION 

Since the first of the Lunar Orbiter missions in August 1966, radar tracking data 
from the Orbiter spacecraft have been intensively analyzed to obtain a better definition 
of the gravitational field of the moon. (See refs. 1 to 6.) These investigations have 
assumed a spherical harmonic expansion of the lunar gravitational potential function simi- 
lar to the one used in this analysis: 

" oo n 

Cq,0+ X Jj (F)X,m ^ (pn,m cos mX + Sn^m sin mx) 
n=2 m=0 _ 

Each of these studies has been an attempt to define a sufficient number of the coef- 
ficients Cn^m Sn^m ^o yield a "good" approximation of the lunar gravitational 

field. The results of these studies have demonstrated that the mathematical representa- 
tion of the observed data is to some extent either erroneous (because of unmodeled 
sources of bias) or incomplete (because of truncation of the potential series). These two 
possibilities are evidenced by systematic residuals in the fit to the tracking data. The 
present analysis has been conducted under the philosophy that the principal inadequacy in 
the mathematical model lies in the representation of localized lunar gravitational anom- 
alies. Hence, the finite -term approximation to the lunar potential function has been 
extended to the thirteenth degree and order in an attempt to provide better resolution of 




these anomalies. Some success has been achieved in the curve -fitting aspects of this 
approach. However, a systematic residual pattern is still evident and points out a 
remaining model incompleteness; possibly, the dominant contribution is from localized 
gravitational anomalies. 


SYMBOLS 

A principal moment of inertia about selenographic X-axis 

a radius of moon, 1738.09 km 

B principal moment of inertia about selenographic Y-axis 

C principal moment of inertia about selenographic Z-axis 

Cn,mjSn,m coefficients of lunar gravitational potential harmonics 
GM product of gravitational constant G and lunar mass M, 4902.58 km^/sec^ 

M lunar mass 

m associated Legendre polynomials of degree n and order m 
r selenographic radius , km 

U gravitational potential function, kmVsec^ 

X,Y,Z selenographic axes 
X selenographic longitude, deg 

<p selenographic latitude, deg 

DISCUSSION AND RESULTS 
Method of Analysis and Data Coverage 

Table I presents approximate orbital characteristics for those phases of the Lunar 
Orbiter missions from which data have been analyzed in a weighted-least-squares , 
differential-correction process. Fundamental to this process is the numerical integration 



of the spacecraft trajectory (Cowell formulation) and the associated variational equations 
utilizing a twelfth-order predictor-correction technique with 13 -digit precision. Initial 
values specified for the spacecraft state and the lunar gravitational coefficients are then 
differentially corrected and the procedure is iterated to achieve a least-squares fit to 
observed data. These observations consist of counting, over 1 -minute intervals, the 
number of cycles in the Doppler frequency shift in a radio signal sent from an earth sta- 
tion to the spacecraft, where it is retransmitted to an earth receiving station. 

The lower part of table I gives a breakdown of the data into the nuiriber of data arcs 
(data representing a continuous spacecraft trajectory), number of data points, and total 
length of coverage per orbiter provided by these arcs. 

The total data currently analyzed amounts to 20 148 coherent two-way Doppler 
observations representing approximately 80 days of tracking. This direct observational 
coverage is illustrated in figures 1 and 2 for the total data and for data taken when the 
spacecraft was below 200 km in altitude, respectively. Because of the line -of -sight char- 
acteristic of the radar tracking system, direct observation of the spacecraft is inter- 
rupted as it moves behind the moon. The selenographic longitude of the spacecraft at 
this interruption depends upon the altitude of the spacecraft above the lunar surface. It 
follows that direct observations relating to the lunar far- side environment must be at 
relatively high altitude, as compared with near-side observations. As a result of radial 
attenuation of gravitational effects, these data are less sensitive to localized gravitational 
anomalies than the near-side data which represent a much wider range of altitudes. 
Comparison of figures 1(b) and 2(b) shows the impact of this altitude constraint on the 
data coverage. Comparison of figures 1(a) and 1(b) with figures 2(a) and 2(b) gives an 
implication of the relative sensitivity to localized gravitational anomalies of near -side 
data compared with far-side data. This high- altitude constraint also applies to a lesser 
degree to the polar regions because of mission design which required periapse photo- 
graphy of near-equatorial sites. 

An additional point with regard to figure 1 is that although areas are shown in which 
no direct observations were taken, these regions are indirectly represented in the data. 
This is true in the sense that the path of the spacecraft over these areas is included in 
at least one of the continuous arcs used in the analysis. Thus, the net effect of the gravi- 
tational forces in these regions is reflected in the spacecraft trajectory and hence in 
the data. However, the absence of direct data and/or low-altitude data remains as an 
important limitation in view of the current incompleteness of the mathematical model. 
Consequences of this limitation are pointed out later in relation to gravitational contour 
maps. 
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Estimates of the Lunar Gravitational Coefficients 

Estimates of the gravitational harmonic coefficients through degree and order 
thirteen are given in table II. The estimate given for Cq q implies a value of 
4902.84 km^/ sec^ for GM. In view of the preceding discussions relating to model 
incompleteness and data coverage, these numbers should be considered as interim esti- 
mates, subject to later refinement. Caution should be exercised in interpreting the num- 
bers in table H, since their validity as gravitational field estimates requires use of the 
complete set of coefficients. For example, based on a sequence of unpublished gravita- 
tional field estimates of different degrees (generated by the first author in analysis of 
varying combinations of data), the field estimates are relatively stable, whereas the dis- 
persion in individual coefficient estimates indicates that they are accurate only to within 
an order of magnitude. Exceptions to this rule are GM and the coefficients C2 Qj ^ 22 ’ 
and C31. Table HI presents values of these parameters. The parameter uncertainties 
listed with the current estimates in column 5 of table HI are based on the dispersion of 
this sequence of values about the current estimates. 

Implications of the Gravitational Field Estimates 

Contour maps of the radial component of acceleration corresponding to the current 
estimate of the lunar gravitational field are given in figure 3. These contours are in 
milligals (10“^ cm/sec2) and represent the deviation in the radial acceleration from that 
of a spherically symmetric field, for which GM = 4902.84 kmS/sec2, at an altitude of 
100 km. The three large Northern Hemisphere positive anomalies occurring below a 
latitude of 40° in figure 3(a) correspond closely to the lunar maria Imbrium, Serenitatis, 
and Crisium. (Refer to the lunar map given in fig. 4.) These three anomalies are in 
substantial agreement with results given in references 7 and 8. Comparison of figure 3(a) 
with figure 1(a) shows relatively dense direct data coverage of these areas with the excep- 
tion of the northernmost portion of the Serenitatis anomaly. This fact, coupled with the 
relative stability of these anomalies throughout the analysis, leads the authors to the con- 
clusion that these contour anomalies are reasonable representations of true lunar gravi- 
tational anomalies. In contrast to these areas, consider the large positive anomaly at 
longitude , -60°; latitude, -35°. Inspection of figure 1(a) shows this anomaly to lie in a 
region devoid of direct observational data. This result may be more than an interesting 
coincidence because if the general "roughness" of figure 3(a) is compared with the general 
data coverage, implied by the combination of figures 1(a) and 2(a), a correlation between 
contour "roughness" and the relative sensitivity of the data to localized gravitational 
anomalies is apparent. If this type of comparison is extended to the lunar far side, the 
correlation becomes even more apparent from the alternating high-low nature of the con- 
tours and their increase in amplitude up to approximately 1000 milligals (1 cm/sec2) in 
the region of 180° longitude. It seems likely that the general "roughness" of the 



gravitational field estimate, although induced by localized gravitational anomalies, is 
more representative of the currently incomplete mathematical model of the data and the 
curve-fitting philosophy of the least-squares method. It is believed that these factors 
lead to a biased estimate in the areas of least sensitive data coverage in order to achieve 
a more valid estimate in the areas of most sensitive coverage. For this reason, no 
attempt has been made to correlate the current estimate of gravitational anomalies in the 
polar and far-side regions with lunar surface features. 

In spite of the low confidence level assigned to most of the individual gravitational 
coefficients, the stability of certain low-degree coefficients permits some reasonable 
inferences to be drawn. In particular, the second-degree coefficients can be related to 
the lunar moments of inertia (ref. 9) under the assumption that (C - A)/B = 0.000629. 
Table IV lists the associated selenographic inertia tensor and the principal moments of 
inertia based on the values in table n. Uncertainties given for the principal moments are 
estimated from the dispersion resulting from use of the sets of second-degree coefficients 
associated with the gravitational field estimates of table HI. These estimates indicate 
that the lunar principal axes differ from the selenographic axes by no more than 5.5° and 
imply an essentially homogenous moon. (C/Ma^ = 0.4 for a homogeneous sphere.) 

The mare -associated gravitational anomalies discussed earlier present a contrast 
to the basic homogeneity implied by the moment-of-inertia estimates. To achieve a 
rough estimate of a possible density distribution causing these anomalies, point-source 
approximations were used for the Imbrium and Serenitatis anomalies. The gravitational 
field of these point sources was expanded in spherical harmonics to the same degree and 
order as the current estimate of the lunar gravitational field. By the method of least 
squares, this expansion was then fitted to the lunar gravitational-field estimate over 
regions in space varying from 50 to 150 kilometers in altitude and lO'^ in selenocentric 
angle from the "axis" of the anomaly. The resulting solutions for the masses and loca- 
tions of the point sources are given in table V and indicate, as one explanation, large mass 
concentrations at depths of approximately 200 and 300 kilometers. Figure 5 illustrates 
the degree of agreement between the point-source solution and the current spherical har- 
monic estimate of the Serenitatis anomaly. Although the gravitational-field analysis is 
not sufficiently stable to justify an exhaustive interpretative analysis, the preceding 
exercise illustrates the type of investigations that can be conducted with a more refined 
gravitational -field estimate and more elaborate source models. 

CONCLUDING REMARKS 

The preceding results represent the current status of a continuing analysis. Prob- 
lem areas associated with mathematical modeling, observational data coverage, and the 
effect of these areas on the validity of the analysis have been alluded too. Specifically, a 
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strong correlation has been shown to exist between the type of data coverage and the 
general "roughness" of the current gravitational-field estimate. This correlation, in 
connection with the instability of the lunar far-side gravitational contours , indicates that 
current estimates of localized gravitational anomalies on the lunar far side are not 
physically valid. However, the stability of the estimates of the near-side gravitational 
contours permits a reliable association of gravitational anomalies with certain lunar 
maria. Also, the stability of the low degree and order harmonic coefficients yields a 
reliable implication that the moments of inertia of the moon are very nearly those of a 
spherical, homogeneous body. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va, April 30, 1971, 
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TABT.-R II.- COEFFICIENTS OF THE THIRTEEKTH DEGREE AM) ORDER SOLUTION 


FOR LimiR GRAVITATIOWAL FIEED 


n 

m 







9 

7 

- 8. 0690242529 

X 

10-11 

8.0754026764 

X 

10-21 


8 

- 1.2649479578 

X 

10-12 

2.6681227345 

X 

10-11 


9 

1.7805101055 

X 

10-12 

- 5.4376209997 

X 

10-13 

10 

0 

- 1.9534030711 

X 

10-5 

0 




1 

3.5640472538 

X 

10-5 

44 . 51233094077 X 

10 "^ 


2 

4.9196310871 

X 

10-9 

5.7169685637 

X 

10-7 


3 

- 2.2358784038 

X 

10 "^ 

- 3.6437585575 

X 

10-8 


4 

3.9141653049 

X 

10-9 

1.1269567323 

X 

10-8 


5 

3.8087666364 

X 

10-10 

- 2.1308938653 

X 

10-9 


6 

- 4.6255135189 

X 

10-11 

1.5652458548 

X 

10-10 


7 

2.7614962503 

X 

10-11 

- 1.5540376682 

X 

10-11 


8 

5.4896743830 

X 

10-12 

- 4.0505042040 

X 

10-12 


9 

9.2407535252 

X 

10-14 

- 3.6136326596 

X 

10-13 


10 

- 7.5223488683 

X 

10-14 

3.3135927064 

X 

10-14 

11 

0 

5.1037452530 

X 

10”5 

0 




1 

- 7.7291781814 

X 

10-6 

2.8479446434 

X 

10-5 


2 

- 3.2273658086 

X 

10-7 

7.2497774428 

X 

10-7 


3 

2.3360853505 

X 

10-8 

5.5966682329 

X 

10-8 


4 

1.0287849331 

X 

10-9 

1.3205804852 

X 

10-8 


5 

- 3.1877004821 

X 

10-10 

- 8.8603014695 

X 

10-10 


6 

6.2535668259 

X 

10-11 

1.2557575513 

X 

10-10 


7 

- 1.3948668302 

X 

10-12 

- 1.6191981587 

X 

10-11 


8 

- 9.6800165098 

X 

10-13 

2.8512390282 

X 

10-13 


9 

- 2.8524797963 

X 

10-13 

2.5679558039 

X 

10-13 


10 

- 1.4876779292 

X 

10-14 

2.0466569257 

X 

10-14 


11 

7.2963987090 

X 

10-16 

- 1.1859856626 

X 

10-15 

12 

0 

9.7617164640 

X 

10-8 

0 




1 

1.2778274406 

X 

10“5 

1.4438415728 

X 

10-6 I 


2 

- 1.5405505424 

X 

10-7 

5.3551638732 

X 

10-7 


3 

5.4002051640 

X 

10-9 

8.7868585746 

X 

10-9 


4 

4.0078293255 

X 

10-9 

1.2737288951 

X 

10-9 1 


5 

9.5030429269 

X 

10-11 

- 4.8376477908 

X 

10-10 1 


6 

1.7093579633 

X 

10-11 

6.0134093338 

X 

10-12 1 


7 

- 7.2263198820 

X 

10-13 

- 5.5677059494 

X 

10-12 1 


8 

1.2421199266 

X 

10-13 

9.7510387449 

X 

10-13 


9 

- 3.7285776082 

X 

10-14 

- 2.4368928834 

X 

10-14 


10 

9.9312843321 

X 

10-15 

- 1.0618795919 

X 

10-14 


11 

5.2033681565 

X 

10-16 

- 4.2870250967 

X 

10-16 


12 

- 2.3855459682 

X 

10-18 

4.5543340194 

X 

10-18 

13 

0 

7.5782274244 

X 

10-8 

0 




1 

- 5.7637379713 

X 

10-8 

2.3610613996 

X 

10-5 


2 

- 1.5991666946 

X 

10-7 

3.4297075982 

X 

10-7 


3 

2.7190638195 

X 

10-8 

9.6757608497 

X 

10-9 1 


4 

7.6732684554 

X 

10-10 

3.0028913747 

X 

10-9 1 


5 

- 2.1117612255 

X 

10-10 

- 1.1279989530 

X 

10-10 


6 

1.8106651096 

X 

10-11 

1.9737118563 

X 

10-11 


7 

- 1.6221201052 

X 

10-12 

- 1.1582995210 

X 

10-12 


8 

- 8.7167027819 

X 

10-14 

8.6768075576 

X 

10-14 


9 

1.0705140705 

X 

10-14 

- 2.7869365761 

X 

10-14 


10 

- 7.5223488683 

X 

10-14 

3.3135927064 

X 

10-14 


11 

- 2,8387572286 

X 

10-16 

2.0316367811 

X 

10-16 


12 

- 9.0639928622 

X 

10-18 

- 1.2565366439 

X 

10-17 


13 

1.1380844737 

X 

10-18 

- 2.6006569878 

X 

10-19 


n 

m 

^ n,m 



Sn,m 



0 

0 

1.00005388274 


— 

- 


2 

0 

- 2.0378761820 

X 

10-4 

0 




1 

1.1051542512 

X 

10-5 

1.3006148287 

X 

10-5 


2 

2.4845211490 

X 

10-5 

- 1.0442661244 

X 

10-8 

3 

0 

2.8439912986 

X 

10-5 

0 




1 

2.4152698939 

X 

10-5 

2.0805931299 

X 

10"5 


2 

7.6322652451 

X 

10-5 

2.2710804784 

X 

10-6 


3 

1.4111967268 

X 

10-5 

- 3.1126414977 

X 

10-7 

4 

0 

3.4688330254 

X 

10-5 

0 




1 

- 1.9891933273 

X 

10-5 

- 8.5100449545 

X 

10-6 


2 

- 2.5443614074 

X 

10-5 

- 4.1656170604 

X 

10-6 


3 

- 5.6321036284 

X 

10-7 

- 2.7682641529 

X 

10-7 


4 

- 5.3981736976 

X 

10-8 

1.1411256460 

X 

10-7 

5 

0 

- 2.6623483022 

X 

10-5 

0 




1 

- 7.4939720619 

X 

10-5 

- 9.2535826576 

X 

10-6 


2 

3.4603507569 

X 

10-7 

3.8762399793 

X 

10-7 


3 

4.0040015934 

X 

10-7 

8.8773988822 

X 

10-7 


4 

9.3924925003 

X 

10-8 

- 1.4794482739 

X 

10-7 


5 

1.8771823966 

X 

10-8 

- 2.1001366939 

X 

10-8 

6 

0 

4.9673195626 

X 

10-5 

0 




1 

- 1.8543318257 

X 

10-5 

- 9.5673010658 

X 

10-6 


2 

- 5.5091031933 

X 

10-7 

1.8138036159 

X 

10-7 


3 

1.1678358571 

X 

10-7 

- 1.1841174024 

X 

10“7 


4 

- 3.3956313198 

X 

10-8 

- 9.7021389127 

X 

10-8 


5 

- 9.6847691920 

X 

10-9 

8.1237744599 

X 

10-9 


6 

- 4.0790023093 

X 

10"9 

2.9739566973 

X 

10-9 

7 

0 

- 7.2781850970 

X 

10-5 

■0 




1 

4.1810263527 

X 

10-5 

- 2.8408094790 

X 

10-5 


2 

- 8.3378247684 

X 

10-7 

- 7.4166186129 

X 

10-7 


3 

- 1.2737535124 

X 

10-7 

- 1.7509754956 

X 

10-7 


4 

- 4.1682647653 

X 

10-8 

- 2.5407115475 

X 

10-9 


5 

1.8658206027 

X 

10-9 

1.0910340047 

X 

10-8 


6 

4.6261404574 

X 

10-10 

- 7.1046094393 

X 

10-10 


7 

3.4334519153 

X 

10-10 

- 1.2401781124 

X 

10-10 

8 

0 

- 1.7411742627 

X 

10-5 

0 




1 

1.3701858067 

X 

10-5 

3.6716642133 

X 

10-6 


2 

4.0795796661 

X 

10-7 

9.1385019071 

X 

10-7 


3 

9.3910161964 

X 

10-9 

- 1.9359154697 

X 

10-7 


4 

- 6.4523074373 

X 

10-9 

2.5059489581 

X 

10-8 


5 

4.5024517927 

X 

10-9 

- 9.1723952804 

X 

10-10 


6 

3.6434678341 

X 

10-10 

- 1.0430917747 

X 

10-9 


7 

- 1.4858612477 

X 

10-11 

2.1525564980 

X 

10-11 


8 

- 3.6012675699 

X 

10-11 

1.1935027479 

X 

10-11 

9 

0 

- 6.5011326253 

X 

10-5 

0 




1 

5.3770904532 

X 

10-7 

- 5.2273480368 

X 

10-6 


2 

- 4.4206541496 

X 

10-7 

- 3.1727246486 

X 

10-7 


3 

- 1.2910912315 

X 

10-7 

- 5.4900413900 

X 

10-8 


4 

- 8.3591988861 

X 

10-9 

2.1240873880 

X 

10-8 


5 

8.7259392768 

X 

10-10 

- 2.5124288158 

X 

10-9 


6 

- 5.1505398116 

X 

10-10 

1.2516041989 

X 

10-10 


9 
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(a) Lunar near side. 
Figure 1.- Total data coverage. 
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(b) Lunar far side. 
Figure 1.- Concluded. 






(a) Lunar near side. 

Figure 2.- Observations taken when the spacecraft was below 200 km altitude. 
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(b) Lunar far side. 


Figure 2.- Concluded. 
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(a) Lunar near side. 


Figure 3.- Estimate of variations in gravitational acceleration (milligals) . 
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(b) Lunar far side. 
Figure 3.- Concluded. 
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Figure 4.- Map of lunar near side. 
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